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An analytical representation for the spatial and temporal dynamics of the simplest of the diffusions 
- Brownian diffusion in an homogeneous slab geometry, with radial symmetry - is presented. This 
representation is useful since it describes the time-resolved (as well as stationary) radial profiles, for 
point-like external excitation, which are more important in practical experimental situations than 
the case of plane-wave external excitation. The analytical representation can be used, under linear 
system response conditions, to obtain the full dynamics for any spatial and temporal profiles of 
initial perturbation of the system. Its main value is the quantitative accounting of absorption in the 
spatial distributions. This can contribute to obtain unambiguous conclusions in reports of Anderson 
localization of classical waves in three dimensions. 
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I. INTRODUCTION 

Diffusion, the propagation of a perturbation in space 
and time, is ubiquitous. It is so ever present that the 
qualitative concept is used even outside the context in 
which it was originally defined, like in the case of spectral 
diffusion or the diffusion of alleles in a population in pop¬ 
ulation genetics. The contemporary approach is to con¬ 
sider diffusion, even in the strict sense of spatio-temporal 
spreadening of an initial perturbation, in a broader con¬ 
text: one then have the classical Brownian diffusion, the 
first diffusion to have a quantitative model for its de¬ 
scription, as a particular case of a whole possible menu 
for diffusion processes. Nevertheless, classical diffusion 
remains the canonical simple model against which all 
other diffusions are gauged against. As a result, diffu¬ 
sions can be labelled as diffusion (per si, the Brownian 
case) or anomalous diffusion [T]. The anomalous case, 
rather than being exceptional, provides the current state- 
of-the-art framework to approach diffusion processes. It 
is subdivided into subdiffusion (slower than the standard) 
or superdiffusion (faster). One of the alternative quanti¬ 
tative descriptions of anomalous diffusion uses fractional 
order derivatives and it is therefore a generalization of the 
classical diffusion equation, although one should state in¬ 
stead that the classical diffusion equation is the special 
case or integer order derivatives BE]- However, no mat¬ 
ter our current understanding of the plethora of diffusion 
processes, the classical, Brownian-Einstein case remains 
our starting point. 

If one uses a stochastic description of diffusion pro¬ 
cesses, one captures the dynamics in the spatial and tem¬ 
poral probability distribution functions, the ones describ¬ 
ing the conditional evolution of the system, given past 
conditions. If the second moment of the spatial distri¬ 
bution and the first moment of its temporal counterpart 
are finite, one can use the classical diffusion equation [2]. 
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This will give a solution, valid asymptotically. Every mi¬ 
croscopic detail can be encapsulated in a meso or macro¬ 
scopic description using an effective diffusion coefficient. 
This coefficient will describe the ensemble, in time and 
space scales large enough when compared to each individ¬ 
ual microscopic interaction. This is the physical insight 
which is the translation of the more abstract concept of 
the central limit theorem. In the classical diffusion case, 
one can encapsulate the microscopies into the mean free 
path, which defines a typical scale for the interaction. 
The physical insight continues to build upon this: it is 
both a typical scale as well as a distance which can be 
used to substitute the whole range of distances for this 
single one. It is this insight which is used below to substi¬ 
tute the whole distances for the initial perturbation, for 
a single individual distance corresponding to the mean 
free path. There is life after a mean free path, though. In 
multiple scattering conditions, after tens to thousands of 
mean free paths, one aims to obtain the description of the 
ensemble, a cumulative response due to a large number 
of (possibly uncorrelated) individual events. It is in this 
context that the diffusion equation is most useful and it 
is for this well developed multiple scattering regime that 
the contribution of this manuscript is adequate. 

The transport of waves is a multiple scattering regime 
is both well known and very important, from the point 
of view of applications. The simplest of the diffusions in 
the title refers to both the simplest case of diffusion-like 
type of interactions, the one amenable to description by 
a diffusion equation, as well as its simplest realization in 
an actual experimental system: a homogeneous slab fi¬ 
nite system. In this system the diffusion is isotropic. The 
historical more important case of this system is perhaps 
the well-known case of plane-parallel stratified stellar at¬ 
mospheres [3]. There is a natural stratification in the 
perpendicular coordinate. It is therefore natural to use 
this coordinate as the ONLY spatial relevant coordinate 
of the system. Eor an actual laboratory setup, this means 
that a plane-parallel incident wave is needed, in order not 
to break the symmetry. Of course, the driving force for 
this is that the full spatial dynamics is amenable to a uni- 
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dimensional formulation [4]. This is useful, to grasp the 
essence, but is of difficult practical realization in a lab. 
The other simplifying assumption which is common is to 
focus only in the stationary state description. A recent 
contribution was made in which the radial profile was 
obtained, in steady state conditions, but not restricted 
to a pseudo-unidimentional formulation [5]. In this con¬ 
text, a full solution of the diffusion equation, describing 
the space and time dynamics, in multiple scattering and 
in space resolved coordinates in axial symmetric condi¬ 
tions (pointlike initial perturbation and isotropic scat¬ 
tering in an homogeneous slab), was attempted. It is an 
analytical series representation for that solution that this 
contribution presents. Under linear system theory condi¬ 
tions, this solution can be used to obtain all the ensemble 
responses to any spatial or temporal initial perturbation, 
by making the respective convolutions. 

As a summary, the simplest diffusion in the ti¬ 
tle refers to: (1) multiple scattering regime, amenable 
to a description by the canonical diffusion equation, 
(2) isotropic interaction in an homogeneous plane- 
parallel slab, maintaining axial symmetry, (3) linear re¬ 
sponse, and (4) scalar waves. 

The manuscript is organized as follows: Secs. II and III 
give an analytic representation for the spatial and tempo¬ 
ral dynamics, first in reflection and transmission (II) and 
after inside the slab (III). Sec. IV presents conclusions. 


for time and space will be s and k± (in this last case, 
conjugate to the radial distance only). 

The diffusion equation describing the propagation of 
the intensity distribution is: 

where the symmetry is emphasized by the explicit 
grouping of spatial coordinates. / is the intensity dis¬ 
tribution, and Dq, a and S the diffusion coefficient, ab¬ 
sorption coefficient and source term. The operator 
refers only to the coordinates defining the perpendicular 
to the symmetry axis. 

For easy reference, / = and = 

jLF, 2 D z^s). I (r_L, 2 :, t) signals axial symmetry. 
jLF, 2 D Laplace-Fourier transform, with 2D em¬ 

phasizing that a 2D spatial Fourier transform is implied. 

The source term is now written S{r±^z,t) = 

S (x) S (y) S {z — (/* -h ^e)) ^ {t) (unit production plus the 
diffusion approximation of the whole of the external ini¬ 
tial perturbation being substituted for a single delta, a 
transport mean free path away from the sample enclosure 
wall; all observables will come out normalized relative to 
the overall strength of the external perturbation). 

In Fourier-Laplace one is left with [ 7 ]: 


II. REFLECTION AND TRANSMISSION 
DYNAMICS 


Fig. [T] shows the notation used in the derivations. 
The 2 :-coordinate defines the coordinate inside slab. /* 
is the transport mean free path, L the slab size, and 
2^6 a so-called extrapolation length (this will be used in 
the boundary conditions for the differential equation; its 
simplest value is given by 2:e = 2/3/*, although more 
elaborate approximations are possible, depending on re¬ 
fraction in the interfaces Hi). Le is an effective width. 
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FIG. 1: Notation used in derivations, z, L, Ze, /*, and Le 
are the slab coordinate and size, an extrapolation length, the 
transport mean free path, and an effective width, respectively. 
Actual slab between coordinates Ze, and Le — Ze- 

For axial symmetric conditions, one can group carte¬ 
sian coordinates into 2 : and r^, the perpendicular dis¬ 
tance to the symmetry axis. The conjugate coordinates 


dz‘^ 
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Using standard Dirichelet boundary con¬ 
ditions 0 m and the extrapolation length, 
{/ {r±,z = 0,t) = 0, / {r±,z = Le,t) = 0} [4], one 

obtains the solution: 


111 sinh \/A (/* + 2 :e) sinh \/A (Te — 2 :) 

TT Dq sinh ^/ALe 

ifz< r ^Ze, (3) 

11 1 sinh (Le— (/*+^e)) sinh v^ 2 : 

TT Do y/A sinh y/AL^ 

if 2 : > /* -h (4) 

with A = k\^^-^. 

The goal is to obtain quantities directly amenable to 
measurement, in simple setups. The sought experimental 
observables are therefore the fluxes (4>) in the exit planes, 
imaged in reflection {R) or transmission (T): 
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T = ^ {z = Le - Ze,t) . 

( 7 ) 


One is left with 
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1 smh Ze) cosh ^/A{Le — Ze) 


= —- 


sinh y/ALe 


. ( 8 ) 


_ 1 sinh v/]4(Le-(/*+ 2 :e)) cosh v/]4(Le-^eLx 
^ ^ sinhViLe 


whose Fourier-Lapace inversion gives: 


1 1 p-at ^ 

n=l 

7* _1_ -y \ 2 2 -Dn* 

t ~F -2^6 \ —71 TT 2 


2 L2 


xn cos ( I riTT- ^ 


,(10) 


for the transmission time resolved profile. For the re¬ 
flection a similar expression is obtained, the only differ¬ 
ence being that the term ( — 1)’^^^ is absent. PLEASE 
NOTE that in all of the remaining manuscript it is shown 
only the Transmission results. The Reflection counter¬ 
parts are obtained simply by omitting this term. This is 
done in order not to overload the manuscript. 

The last equation is obtained by Laplace inverse trans¬ 
form in the complex plane (Bromwich integral with 
poles at ^/ALe = iriTr) and Eourier inverse transform 
which, for axially symmetric functions, takes the simple 
form (r) = (k) Jq (kr) kdk^ using the zero 

order Bessel function [iQiin]. 

Three further important observables are to¬ 
tal intensity, integrated in space, j^total _ 

T {r A) r dO dr ^ the steady state radial profiles, 
j^ss rjn total signal also in 

steady state, obtained by spatial integrating the radial 
profiles 7 ^‘S'S',total _ j^ss ^ dOdr: 


TTOTAL 

n=l 

^e\ . f F Zq\ —n^TT^ \ 

xn COS ( nTT— I sm | nn —^- | e , (11) 
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Q +00 
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^ (^ + Don+ ) 
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where Kq is the modified Bessel function of the second 
kind. 


Eqs. ( |T0p3l ) constitute the main contribution of this 
manuscript. Compare also the time resolved total (not 
spatial resolved) fluxes in refs. [12] and [T3] . 

Please note that, without absorption {a = 0) 


J^SS,TOTAL ^ jiSS,TOTAL ^ 

= - + --- cos ( (2n + 1) TT— ) X 

7r+'2n+l V LeJ 


77 = 0 


X sin ( (2n + 1) tt— ^ ^ ) = 1, 


(14) 


which express energy conservation (and is therefore a 
simple consistency check). 

Eqs. ([T^) are best recast in dimensionless coordi¬ 
nates. The characteristic scales are: 


L’ 

(15) 
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and the previous results become: 
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n cos J 

Ko (+ n^TT^ + a'r') , 


( 20 ) 


+ 00 
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Two aspects are worth mentioning. The first is the 
asymptotic behaviour, in space and time. This is ob¬ 
tained in the previous equations making n = 1, which is 
the slowest decaying term. The asymptotics are: 
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(22) 

t' 

rj^TOTAL 

oc 

(23) 

tSS 

OC Aq . 

(24) 


The time-resolved distributions are Gaussian (with 
variance linear with time). The total signal, integrated 
in space, is exponential with t'. And the steady state 
radial distributions are quasi-exponential with the radial 
distance r'. This last statement is due to the facts mm- 
(1) Ka (x) = (ix), for —tt < arg{x) < tt, 

and (2) Kq + a'r') oc r'~^^‘^exp (—+ a'r'). 
The has a much weaker dependence on the dis¬ 

tance than the exponential, and a approximate exponen¬ 
tial dependence follows. 

The second aspect is the mean displacement, as can 
be judged from the radial profiles in reflection or trans¬ 
mission. It is well known that Brownian diffusion gives 
a mean square displacement that is linear with time in 
two dimensions. In fact, the time dependence of this 
mean square displacement is the most important param¬ 
eter distinguishing classical from anomalous diffusion. It 
is therefore tempting to estimate the mean radial dis¬ 
placement^ relative to the symmetry axis, in the imaged 
fluxes. The result is (r'^) oc t', showing that, the physi¬ 
cal signature of classical diffusion survives in the imaged 
fluxes for the three dimensional slab (see also misiiis]). 

Figs. [^to|4 provide a graphical illustration of the re¬ 
sults in eqs. (18[^). Fig. shows the radial profiles, 
in steady state, illustrating the quasi-exponential radial 
asymptotic. This asymptotic regime was recently veri¬ 
fied experimentally 0. Fig. I shows time resolved total 
transmission and reflection, illustrating also the asymp¬ 
totics of eq. (23). Fig. shows that the series in the 
previous equations reproduce the trend known for the 
diffusion in the slab: overall transmission in steady-state 
in inversely proportional to slab width, in the absence 
of absorption [4]. It further shows that this scaling is 
broken, if absorption is present. 

One important aspect in the derivations was that the 
slab was considered infinite, in the transverse direc¬ 
tions. The finite size could destroy the quasi-exponential 
asymptotics of the radial profiles. Reference [5] gives 
empirical evidence that this asymptotic is amenable to 
experimental measurement, for an actual slab represen¬ 
tative of simple experimental conditions (aspect ratio of 
slab diameter to width roughly 10). 

The spatial profiles in fig. [^needed special care in their 
evaluation near the origin. The modified Bessel function 
of the second kind diverge at the origin. Although it 
is not provided an analytical proof of the series conver¬ 
gence, an empirical study of the series convergence was 
nevertheless conducted. Using a sufficient high number 
of terms, it was verified in all cases the series convergence. 
The convergence is however very slow. The fixed number 



FIG. 2: (Color online) Steady state radial profiles, in trans¬ 
mission (T) and reflection (R), using eq. ( [20| (2000 terms, 
r = 1, Le = 50, Ze = 2l3l*). o' is a normalized absorption 
coefficient. 



FIG. 3: (Color online) Time resolved total intensity in trans¬ 
mission (T) and reflection (R), using eq. (20) (2000 terms, 
r = 1, Le = 50, Ze = 2/3/*). a' is a normalized absorption 
coefficient. 


of terms in fig. correspond to a conservative estimate, 
based in the smallest distance considered (see also 0 ). 


III. INSIDE OF SLAB DYNAMICS 


The results presented thus far are focused in the dy¬ 
namics of the transmission and reflection fluxes, since 
these are the ones most easily measured. The dynamics 
inside the slab is nevertheless of importance. The proce¬ 
dures used in the inverse Laplace and Fourier transforms 
of eqs. (|8][^, can be equally applied to obtain the dynam¬ 
ics inside the slab. Eqs. (H give, first for the dynamics 
in direct physical variables 
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FIG. 4: (Color online) Trnasmitance (T), reflectance (R) and 
sum (T+R), as a function of inverse efective slab width, using 
eq. (21) (2000 terms, /* = 1, ; 2 ;e = 2/3/*). a' is a normalized 
absorption coefficient. Without absorption, R + T = 1 which 
is a manifestation of energy conservation. 
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(30) 


Finally, fig. represents the steady-state intensity in¬ 
side the slab, integrated over radial distances, as a func¬ 
tion of spatial coordinate. This result encompasses the 
known trend, for a unidimensional formulation of diffu¬ 
sion in absence of absorption: the intensity is linear in the 
difference to a coordinate inside the slab, situated one 
mean free path away from the boundary (this position 
comes from the approximation in the source term, the 
linear dependence comes from diffusion itself) [121 US]- 
It further quantifies to what extend absorption destroys 
the linear dependence with distance. The inset gives evi¬ 
dence that, at least for strong absorption, the dependence 
with distance can be reasonably approximated by an ex- 
^ponential (see also [4]). 



X 


and, then for its reformulation in dimensionless coor¬ 
dinates 


FIG. 5: (Color online) Steady state intensity distribution 
inside slab, using eq. ( |30t (2000 terms, /* = 1, Le = 50, 
;2;e = 2/3/*). Total intensity, integrated in a plane perpen¬ 
dicular to the z coordinate, a' is a normalized absorption 
coefficient. Actual slab between coordinates and Le — Ze. 
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IV. CONCLUSIONS 

This contribution presents an analytical representa¬ 
tion for the full dynamics in the simplest of the possi¬ 
ble realizations of a diffusion process: classical diffusion 
of scalar waves in a fully developed multiple scattering 
regime, in an homogeneous plane-parallel slab. These 
conditions preserve axial symmetry, which was exploited 
in the derivations. The analytical series is strictly valid 
for a unit initial perturbation. However, under linear sys¬ 
tem response conditions, the convolution with the spatial 
and temporal distributions of the perturbation allow a so¬ 
lution for all other cases. Although this is conceptually 
trivial, even a numerical implementation of a 2D spa¬ 
tial convolution of two radial symmetric functions (initial 
external perturbation and system response) must imple¬ 
ment a proper definition of the convolution as a two- 
dimensional convolution m- For other alternative ap¬ 
proaches for taking into account the finite dimensions of 
the external perturbation see refs. [11 [17]. Ref. [TT] can 


additionally provide a basis for a generalization of this 
contribution to non axial symmetric situations. 

It is expected that this contribution stimulates quan¬ 
titative detailed comparison with experimental results, 
obtained in actual setups not amenable to a realistic uni¬ 
dimensional formulation of the geometry. Its main con¬ 
tribution is the derivation of the spatial resolved profiles, 
including absorption. These can help in obtaining unam¬ 
biguous conclusions in reports of Anderson localization of 
classical waves in three dimensions, a subject of current 
interest and debate [I5l[li[l8]. 
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